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Abstract. The physical origin of high-frequency QPOs (HFQPOs) in black- hole X-ray binaries 
remains an enigma despite many years of detailed observational studies. Although there exists 
a number of models for HFQPOs, many of these are simply "notions" or "concepts" without 
actual calculation derived from fluid or disk physics. Future progress requires a combination of 
numerical simulations and semi-analytic studies to extract physical insights. We review recent 
works on global oscillation modes in black-hole accretion disks, and explain how, with the help 
of general relativist ic effects, the energy stored in the disk differential rotation can be pumped 
into global spiral density modes in the disk, making these modes grow to large amplitudes 
under certain conditions ( "corotational instability"). These modes are robust in the presence 
of disk magnetic fields and turbulence. The computed oscillation mode frequencies are largely 
consistent with the observed values for HFQPOs in BH X-ray binaries. The approximate 2:3 
frequency ratio is also expected from this model. The connection of HFQPOs with other disk 
properties (such as production of episodic jets) is also discussed. 
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1. Introduction 

High-Frequency Quasi-Periodic Oscillations (HFQPOs) in the X-ray fluxes have been 
observed from a number of black-hole (BH) X-ray binaries since the late 1990s. Their 
frequencies (40-450 Hz) are comparable to the orbital frequencies at the Innermost Stable 
Circular Orbit (ISCO) of BHs with mass M ~ lOM©, and do not vary significantly in 
response to sizeable (factors of 3-5) changes in the luminosity. Compared to the low- 
frequency QPOs, the HFQPOs have low amplitudes (0.5 — 2% rms at 2-60 keV) and low 
coherence {Q ^ 2 — 10), and are only observed in the "intermediate state" (or "steep 
power law state") of the X-ray binaries (for reviews, see Remillard & McClintock 2006; 
Belloni et al. 2012). In addition, several systems show pairs of QPOs with frequency ratios 
close to 2 : 3. It is generally recognized that HFQPOs may provide important information 
about the BH (mass and spin) and about the dynamics of inner-most accretion flows. 

Current status of QPO Theory/Models: The physical origin of HFQPOs is cur- 
rently unclear. Below we comment on some of the ideas/models that have been suggested. 

(i) In the spot models (Stella et al. 1999), the HFQPOs arise from the Doppler modu- 
lation of the radiation from spots orbiting in the inner part of the BH accretion disk. It 
is not clear how the position of the spot (a free parameter in these models) is determined 
and how the spot can survive the differential rotation of the disk. 

(ii) Another class of models identify HFQPOs with various oscillation modes of a finite. 
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pressure-supported accretion torus (Rezzolla et al. 2003; Blaes et al. 2006). It is not clear 
that the accretion flow can be approximated by a torus and how the position and size of 
the torus are determined. 

(iii) The harmonic relation between the observed frequencies led Abramowicz & Kluz- 
niak (2001) to suggest that the HFQPOs are a result of a nonlinear resonance. However, 
so far no fluid dynamical model producing these resonances has been developed. Detailed 
calculations of the resonant coupling between the epicyclic modes in slender tori indicates 
that such resonance is very weak (Horak 2008). 

(iv) A large class of theoretical models is based on the relativistic diskoseismology 
(Kato 2001). The ^-modes (also called inertial-gravity modes) have attracted most at- 
tention because their existence does not require a reflective inner disk boundary. These 
oscillations have at least one node in the vertical direction and their restoring force re- 
sults from rotation and gravity. Unfortunately, the non-axisymmetric ^-modes are either 
damped due to corotation resonance (Kato 2003; Li et al. 2003) or have frequencies too 
high compared to HFQPOs. The axisymmetric ^-mode trapped around the maximum 
of the radial epicyclic frequency may account for some single HFQPO frequencies (see 
Wagoner 2012), and several studies suggested that they may be resonantly excited in the 
warped or eccentric disks (Kato 2008; Ferreira & Ogilvie 2008). However, the self-trapping 
property of ^-modes can be easily destroyed by a weak (sub-thermal) disk magnetic field 
(Fu & Lai 2009) and turbulence (Arras et al. 2006; Reynolds & Miller 2009). 

Related to (iv) above, perhaps more promising is the disk p-modes. In the following, 
we review recent works on disk p-modes from the Cornell group and other related efforts 
(see Lai & Tsang 2009, Tsang & Lai 2008,2009a,b,c; Fu & Lai 2009,2011a,b,2012a,b; 
Horak & Lai 2012). 

2. Disk P-modes: Corotational Instability 

Disk j9-modes (also called inertial acoustic modes) represent nearly horizontal oscilla- 
tions with almost no vertical structure. They are trapped between the inner boundary 
of the disk and the inner Lindblad resonance (ILR; see below) Because of their simple 
2D structure, the basic wave properties of the p-modes (e.g., propagation diagram) are 
not strongly affected by disk magnetic fields (Fu & Lai 2009) and are likely robust in 
the presence of disk turbulence. Lai & Tsang (2009) (see also Tsang & Lai 2008,2009c) 
showed that the non-axisymmetric p-modes can naturally grow due to the corotational 
instability. This can be understood as follows. 

Consider a 2D fluid disks without self-gravity and magnetic field. For the disturbances 
of the form e*^^~*^^, where m > and ou is the wave (angular) frequency, the WKB 
dispersion relation for density waves reads 

uj'^ = -\- k'^cl, (where uj = co — mft) (2.1) 

where Q{r) is the disk rotation frequency, hz{r) is the radial epicyclic frequency, kr is the 
radial wavenumber, and Cg is the sound speed. Thus waves can propagate either inside 
the Inner Lindblad resonance (ILR, defined by a; = — ac) or outside the Outer Lindblad 
resonance (OLR, uj = k). The region between ILR and OLR is evanescent ("corotation 
barrier") except for a narrow "Rossby zone" near the corotation resonance (CR, where 
oj = 0). The significance of CR is that waves can be absorbed by the background flow 
since the wave patten speed uj/m matches the flow rotation rate ft at CR. Depending on 
the sign of dC,/dr^ where C, = j (2r]E) is the vortensity of the background flow (E is the 
surface density of the disk), the Rossby zone can lie just inside or outside CR, and the 
corotational wave absorption can be either positive or negative (see Fig. 1). 
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This corotational instability of non-axisymmetric p-modes arises because of two effects: 
(1) Since the waves inside the ILR carry negative energies while those outside the OLR 
carry positive energies, the leakage of the p- waves through the corotation barrier (between 
ILR and OLR) leads to mode growth. (2) More importantly, when the vortensity of the 
disk flow, ( = k!^ / {2QTi) [this applies to barotropic fluid, for which the pressure is a 
unique function of the density; see Tsang & Lai (2009c) for the non-barotropic case], 
has a positive slope at the corotation radius wave absorption at the corotation resonance 
leads to amplification of the trapped p-mode. 

General Relativity plays an important role in the corotational instability of p-modes. 
In Newtonian theory, k = fl^ so d(/dr < when S does not vary strongly and hence the 
mode cannot grow. In GR, reaches a maximum at a finite radius before dropping to 
zero at nsco- Thus, if the mode frequency is such that d(/dr > at corotation, then 
the mode can grow. 

3. Overstable P-modes and HFQPOs 

Our linear calculations and nonlinear simulations (based on pseudo-Newtonian poten- 
tial) have demonstrated the corotational instability of disk p-modes (Lai & Tsang 2009; 
Tsang & Lai 2009c; Fu & Lai 201 la, 2012a). In particular, the low-order p-modes (trapped 
between the inner disk edge and the ILR) have the following properties: 

(i) The mode frequency is /^ml^in, where P = 0.5-0.75 (depending on the disk 
models and the inner boundary condition. The mode frequency varies (by about 10%) 
as M changes by a factor of 3. The frequency is consistent with the estimated BH mass 
and spin of X-ray binaries. The frequency ratio is approximately 1 : 2 : 3, • • •, but not 
exactly. 




Figure 1. Super-reflection and corotational instability of nonaxisymmetric waves in rotating 
disks. The solid curves depict the effective potential for wave propagation: Density waves can 
propagate either inside the ILR or outside the OLR. Between the ILR and OLR lies the coro- 
tation resonance (OR) and the associated Rossby zone (where the effective potential is neg- 
ative). An incident wave e"*^'^^ propagating from small radii towards the corotation barrier 
will be partially reflected (IZe^^'^^) and partially transmitted (Te*^'^^). Since the waves inside 
CR (r < tcr) carry negative energies, while those outside CR carry positive energies, we have 
(-1) = (-l)|7^|^ + |r|^ + Vc from energy conservation, or |7^|^ = 1 + |T|^ + Pc, where Vc is 
the wave energy absorption at CR. When the background disk has dQ/dr > at CR [where 
C = Hi'^/{2QTi) is the vortensity], the Rossby zone lies just outside CR, thus Pc > and 
|7^|^ ^ 1 -\- T>c > 1, we have super-reflection (the left panel). If the disk has a reflective in- 
ner boundary (represented by the red vertical line), then the stationary waves (modes) trapped 
between nn and ILR will be overstable (with growing amplitude) . When the background disk has 
d(^/dr < 0, there is no super-reflection (|7^| < 1) and the modes are damped by the corotational 
wave absorption (the right panel). 
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(ii) The mode growth rate (due to corotation resonance) is much smaher for m = 1 
than for m = 2,3,4, •• •. Since high-m modes are hard to observe, it is reasonable to 
expect that the most visible modes are m = 2,3. 

(iii) The nonlinear mode frequency is close to the linear mode frequency (Fu & Lai 
2012b). 

4. Issues, Complications and Prospects 

All the above make p- modes a promising candidate for explaining HFQPOs in X- 
ray binaries. However, there are a number of issues and complications that must be 
considered: 

(i) Mode Damping due to Radial Infall: For standard thin disks, the radial infall 
velocity increases rapidly near nsco - this can be described in the slim disk model. 
Such infall tends to damp the p-modes (but not completely, because of the sharp density 
gradient around risco; see Lai & Tsang 2009). So we have a competition between mode 
growth (due to corotation resonance) and damping (due to infall). The net result is 
that for standard thin disks, p-modes are likely damped. This is consistent with the 
observations that no HFQPOS are detected in the thermal (low-soft) state of BH X-ray 
binaries. 

(ii) Effects of magnetic fields: However, real disks (particularly those in the "inter- 
mediate state" when HFQPOs are observed) are more complicated, and have magnetic 
fields. As noted above, the p-mode frequencies are only slightly affected by the magnetic 
field, provided that rin is close to risco- However, the mode growth rates can be sig- 
nificant modified by the B-fieid. There are three effects: (1) Toroidal fields inside the 
disk tend to surpress the corotational instability (Fu & Lai 201 la, b). (2) Magnetic fields 
can accumulate inside nsco, forming a magnetosphere around the BH. Thus, the inner 
disk edge may be more reflective than the standard thin disk (Tsang & Lai 2009b; Fu & 
Lai 2012a). This can enhance the net growth rate of p-modes. (3) Large-scale poloidal 
magnetic fields threading the disk can provide a coupling between the disk and the disk 
magnetosphere (Yu & Lai 2013, in prep; see Yu & Lai 2012 for a related problem). This 
tends to enhance the corotational instability and increase the growth rate of p-modes 
(see also Tagger & Varniere 2006). 

The last two points may be of particular importance. Large-scale magnetic fields pro- 
vide a natural means to launch jet s/outfiows from accretion disks. Magnetic fields thread- 
ing the BH can also lead to jets from the rotating BHs (McKinney et al. 2012). In BH 
x-ray binaries, episodic (ballistic) jets are observed in the intermediate state, and this is 
the same state during which HFQPOs are detected. In our model, the disk and corona 
(coupled by a large-scale poloidal 5-field) oscillate together, with the "clock" mainly set 
by the disk dynamics. The corona is important, however, since it can lead to the variation 
of hard x-ray photon flux as observed in HFQPOs. 

The current understanding of accretion flows in the intermediate state is rather poor 
(Done et al. 2007). It is even not clear that rin extends to nsco in this state. We suggest 
that in the intermediate state, large-scale B-fields are temporarily/episodically created 
(e.g., magnetic fields buoyantly rise from the collapsed ADAF state and followed by 
reconnection) . This will then allow the p-modes to grow and manifest as HFQPOs; at 
the same time, it will lead to the production of episodic jets. 

(iii) Full GR Calculations. Horak & Lai (2012) presented full GR calculations of 
overstable p-modes of BH accretion disks. In particular, they derived the GR version 
of the corotational instability criterion and generalized the Newtonian disk vortensity. 
The results are qualitatively in agreement with the pseudo-Newtonian results discussed 
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above, with one important difference: For Kerr BHs with the spin parameter close to 
unity, the mode frequency is close to mQin (i.e., /3 ^ 1). Assuming rin = nsco, this would 
lead to a tension between the simplest p-mode theory and the observational constraints 
on the BH mass and spin. The discrepency becomes most severe for GRS 1915+105 
(M = 14zb4.4M0), which may have two pairs of HFQPOs (41 and 67 Hz, 113 and 168 Hz), 
and whose spin parameter has been constrained to be a > 0.975 using the continuum 
fitting method (see Narayan & McClintock 2012). Thus, it appears that for the simplest 
disk models, non-axisymmetric p- modes have frequencies that are too high compared to 
the observation of HFQPOs in BH X-ray binaries (assuming the BH parameters obtained 
using the continum fitting method is reliable), at least for some systems. Several effects 
may decrease the theoretical p-mode frequencies. For example, a higher disk sound speed 
leads to lower mode frequencies, and a steeper surface density profile (larger p in S oc r~^) 
reduces cJmin, the minimum frequency for the corotational wave absorption to amplify 
the mode. Magnetic fields may also play an important role [see (ii) above]. In particular, 
our calculations so far assume that the inner disk radius coincides with the ISCO. This 
may not be the case during the "intermediate state" (Done et al. 2007; Oda et al. 2010). 
As noted before, a major uncertainty in calculating the disk p-modes is the inner disk 
boundary condition. When magnetic fields advect inwards in the accretion disk and 
accumulate around the BH, (Bisnovatyi-Kogan & Ruzmaikin 1974; Igumenshchev et 
al. 2003; Rothstein & Lovelace 2008; Guilet & Ogilvie 2012), the inner disk radius rin 
may be larger than nsco- This reduces relative to Hisco, leading to lower p-mode 
frequencies. 

(iv) Numerical Simulations. Wen & Lai (2012b) carried out 2D inviscid hydrody- 
namic simulations of overstable p-modes in BH accretion discs assuming that the disk 
inner boundary is sufficiently reffective. The simulations confirmed the linear growth of p- 
modes and showed that the mode growth saturates when the radial velocity perturbation 
becomes comparable to the disc sound speed. During the saturation stage, the primary 
disc oscillation frequency differs only slightly (by less than a few percent) from the linear 
mode frequency. Sharp features in the fluid velocity profiles at this stage suggest that 
the saturation results from nonlinear wave steepening and mode-mode interactions. 

Such 2D simulations, while useful, do not capture various complexities (e.g., magnetic 
field, turbulence, radiation) associated with real BH accretion discs. Numerical MHD 
simulations (including GRMHD) are playing an increasingly important role in unravel- 
ing the nature of BH accretion fiows. A number of "numerical QPOs" have been reported 
in the literature. For example, Henisey et al. (2009,2012) found evidence of excitation of 
wave modes in simulations of tilted BH accretion disks. Hydrodynamic simulations using 
Of- viscosity (Chan 2009; O'Neill et al. 2009) showed wave generation in the inner disk 
region by viscous instability (Kato 2001). The MHD simulations by O'Neih et al. (2011) 
revealed possible LFQPOs due to disk dynamo cycles. Dolence et al. (2012) reported 
transient QPOs in the numerical models of radiatively inefficent fiows for Sgr A"*". McK- 
inney et al. (2012) found QPO signatures associated with the interface between the disc 
inflow and the bulging jet magnetosphere (see Fu & Lai 2012a). It is not clear if these 
QPOs are directly comparable with the observations of BH X-ray binaries. Obviously, 
much work remains to be done to understand the physics of HFQPOs. A combination of 
full numerical simuations and semi- analytic works (together with "controlled" numerical 
experiments), as well as future high signal-to- noise observations of HFQPOs (e.g., with 
LOFT; Feroci et al. 2012), are needed to make progress in this field. 
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